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In previous work we have found a regime in ballistic quantum dots where interelectron interactions 
can be treated asymptotically exactly as the Thouless number g of the dot becomes very large. 
However, this work depends on some assumptions concerning the renormalization group and various 
properties of the dot obeying Random Matrix Theory predictions at scales of the order of the 
Thouless energy. In this work we test the validity of those assumptions by considering a particular 
ballistic dot, the Robnik-Berry billiard, numerically. We find that almost all of our predictions based 
on the earlier work are borne out, with the exception of fluctuations of certain matrix elements of 
interaction operators. We conclude that, at least in the Robnik-Berry billiard, one can trust the 
results of our previous work at a qualitative and semi-quantitative level. 

PACS numbers: 73.50.Jt 



I. INTRODUCTION 

The quest for a satisfactory theory of quantum dots is 
driven not only by their obvious importance as mcso- 
scopic devices revealed by a series of groundbreak- 
ing experiments 0, but also by their challenge as a 
unique confluence of disorder, interactions and finite- 
size effectsH. For weak interactions, the Universal 
Hamiltonian[3, Q (UH) provides a satisfactory descrip- 
tion. For ballistic/chaotic quantum dots, we have 
es pousedH H an approach based on the fermionic 
Renormalization Group (RG), 1/iV expansions and the 
fact that energy eigenstates around the Fermi energy in 
disordered systems ouriit to be described by Random 
Matrix Theory (RMT)0, Our approach not only 

explains the UH as a fixed point of the RG but also de- 
scribes the physics outside its basin of attraction. It pre- 
dicts a phase transition at strong coupling and allows a 
fairly detailed study [71 of the new phase and the quan- 
tum critical region[ll| separating it from that governed 
by the UH. 

Our results, however, were predicated on a variety of 
RMT and RG assumptions. To test our assumptions 
and the conclusions deduced from them, we performed 
a detailed numerical study on a ballistic but chaotic bil- 
liard (the Robnik-Berry billiard [l2|) and we report our 
findings here. Our expectations are borne out, with one 
notable exception. 

We recall our strategy and assumptions briefly so that 
the reader may see in advance what sort of ideas are put 
to test in our study. In the primordial problem of interest 
to us one has electrons confined to a ballistic dot of size 
L, with no impurities inside, and edges so irregular that 
classical motion is chaotic. The electrons experience the 
Coulomb interaction. In momentum space, all momenta 
within the bandwidth (of order kp, the Fermi momen- 
tum) exist. The semiclassical ergodicization time for an 
electron within the dot is a few bounces, or Terg — L/vp. 



By the uncertainty principle this leads to an important 
energy scale, the Thouless energy Ex — Hvp/L which 
has a dual significance. First, it controls the dimension- 
less conductance of a dot strongly coupled to leads, as 
follows. Since the transport through the dot takes place 
in a time such that energy is uncertain by an amount 
Et, all single particle states that fit into this band will 
each contribute a unit of dimensionless conductance. If 
the average single-particle level spacing is S, then the di- 
mensionless conductance is g = Et/S. Second, in the 
other limit of dots very weakly coupled to leads (which 
we focus on in this work), the Thouless band of width 
Et centered on Ep, marks the scale deep within which 
RMT should apply to the energies and eigenfunctionsQ- 
In this context g is better denoted the Thouless number. 

Since we are only interested in a narrow band of ener- 
gies of width Ep << Ep, the first step in the program^ 
is to use the RG for fermionsQ to get an effective low en- 
ergy theory by eliminating all momentum states outside 
Ep- Should we worry that we are not eliminating exact 
single particle eigenstates (labeUed here by a)? No, be- 
cause the disorder due to the boundaries will mix momen- 
tum states at roughly the same energy, and it does not 
matter whether we eliminate momentum states within 
any annulus of energy thickness Ep or the single particle 
states they evolve into. Indeed, even the mixing within 
Ep is due to the fact that momentum itself not well de- 
fined in a finite dot, a point we will elaborate on shortly. 
However, once we come down to within Ep of Ep, we 
cannot eliminate the remaining states in one shot since 
it is the flow of couplings within this band that is all 
important in the RG. 

Now it is knownQ that the clean system RG (justified 
above) leads to Landau's Fermi liquid interaction^^ 



V = J2 ^(^k - 9u')5n{k)5n{k') 



(1) 
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at an energy scale Ep which is small compared to Ep. 
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But since is a bulk scale it can always be made 
larger than Ej which vanishes as L — *■ oo. Thus Murthy 
and MathurQ perform their RG on the hamiltonian (fo- 
cussing on the spinless case for simplicity) 



H 



(2) 



where 
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x(0^(k')05(k)-0,-(k')07(k)) 



(3) 



is simply the Landau interaction written in the basis of 
exact eigenstates, a statement that needs some elabo- 
ration. In usual RMT treatments, 0Q(k) is the exact 
eigenstate a written in the infinite dimensional basis of 
all momentum states. In our version which uses the 
RG to reduce the Hilbert space, the states labeled by k 
are approximate momentum states with an uncertainty 
(5fc ~ 1/i in both directions. The number of such wave 
packets that fit into an annulus of radius hp and thick- 
ness Et/vf is 0{kpL) = g. We call them the Wheel- 
of- Fortune (WOF) states, see Figure (0. One way to 
construct such packets is to pick g plane waves of equally 
spaced momenta on the Fermi circle and to chop them 
off at the edges of the dot to respect the boundary con- 
ditions. This is what we mean by k in ^^(k) = (k|a). 




Assumption II: The energy eigenvalues obey 
RMT statistics as do the wavefunctions. For example 
we assume that the ensemble averages (denoted by ()) 
obey 

(0;(ki)0^(k2)0:(k3)<^.(k4)) = ^ + 0(1/5^) (4) 

To some, Assumption I seems remarkable - How can 
we furnish in advance, independent of the dot shape a 
basis of g states for expanding the g exact eigenstates 
within After all, these eigenstates are supposed to 
resemble those of a random matrix. The point is that no 
matter how chaotic the dot, it can only mix states at the 
same energy. While this sounds like Berry's ansatz^^l it 
is somewhat different in both scope and content: Berry's 
ansatz states that every exact eigenstate a can be ex- 
panded in terms of an infinite number of k states (with 
the same energy ea = k^/2m)in the bulk, while we claim 
that only g of WOF states are needed. Secondly, we claim 
that the same g WOF states can be used to describe all 
the states within the Thoulcss band. 

In this work we will show that these states are nearly 
orthonormal and that the exact state right in the middle 
of the band has 99.9% overlap with the WOF states. The 
success of this extension of Berry's conjecture to a finite 
dot exceeds our expectations in this regard. However, 
we find that the g WOF states become less effective at 
describing states as we move away from Ep- the overlap 
drops to 50% at E^ = Ep ± Et/2. This prepares us 
for the possibility that nonuniversal quantities may be 
quantitatively inaccurate in our approach. 

As for Assumption II, we have verified RMT behav- 
ior for the eigenvalues (as have others before us [Til H^ 'l 
but not the eigcnfunctions. What we did instead was 
to see what extent the solution of a specific dot resem- 
bled the picture we drew based on these two assump- 
tions. We begin by describing how one starts from Eqn. 
0, which describes the effective hamiltonian, and use 
our two assumptions with large-A^ ideas to make our 
predictions!^ Q. These predictions are asymptotically 
exact as g ^ oo. 

First one expands the Landau function as 
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FIG. 1: The g wheel-of-fortune states in the Thouless band. 
They are packets in k space with average momenta equally 
spaced on the Fermi circle. 

We are now ready to state our two assumptions. 
Assumption 1 ; We assume that the g approximate mo- 
mentum states generated as above form a complete basis 
for the g exact eigenstates in the Thouless shell. 



(5) 



Barring accidents, the phase transition occurs in one 
channel with some m (recall superconductivity). This 
allows us to focus on a single Um ^ 0, ignoring all oth- 
ers. Then we carry out a Hubbard-Stratovich transfor- 
mation on the interaction using a collective field cr. We 
then formally integrate out the fermions and get an ef- 
fective action S{(t) for cr. In this process we make use 
of assumption II. The action in terms of cr is obtained 
by summing one loop Feynman diagrams with varying 
numbers of external cr's connected to a single fermion 
line running around the loop. Each diagram is a sum 
over fermion energy denominators multiplying products 
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of a string of 0Q.(k)'s. We are able to show that these 
products may be replaced by their ensemble averages in 
the large g limit. In other words the sum over so many 
terms in each diagrams leads to self-averaging. For the 
averages we use relations like Eqn. When this is 

done, the effective action can be cast into a form which 
has a in front of itj^Q, so that the saddle point gives 
exact answers as g — > oo. 

At this point let us collect all the results and predic- 
tions of the RMT + large-iV theoryQ with a view to 
comparing them with similar results without using As- 
sumptions I and II on the Robnik-Bcrry billiard. 

• In the large-g limit there is a sharp transition to 
a phase in which cr acquires a vacuum expectation 
value. The critical value of Um in our approxima- 
tion turns out to be — 1 / log 2 in the spinlcss case 
and —1/2 log 2 in the spinful case. The true critical 
value is most likely to be the bulk value —2 (spin- 
less) or — 1 (spinful) , as has been found in an explic- 
itly solvable model by Adam, Brouwer, and Sharma 
recently This is an example of the nonuniversal 
quantity alluded to earlier, that we cannot predict 
exactly in our approach even as <? ^ oo. 

• For finite 5, instead of a sharp phase transition, 
there is a crossover from the weak-coupling regime 
through a quantum critical regime to a strong- 
coupling regime. Due to the explicit symmetry- 
breaking at order 1/17, there is always some nonzero 
order parameter, which increases to a number of or- 
der g (in the normalization we use here, which is 
different from that of ref. ) in the strong-coupling 
regime. 

• For symmetry-breaking in odd angular momentum 
channels there are two exactly degenerate minima 
for every sample arising from time-reversal invari- 
ance. 

• The ground state energy at the minimum in the 
strong-coupling regime is lower than that in the 
weak-coupling regime by a number of order g^5. 

• The effective potential landscape in the strong- 
coupling regime is in the approximate shape of a 
Mexican Hat, with the ripples at the bottom of the 
hat being of order g5. 

• In the quantum-critical and strong-coupling 
regimes, even low-energy quasiparticles acquire 
large widths given on average by 



r(£)«-log(e/5) 

TT 



(6) 



We found that most of these predictions are verified 
by our numerical results on the Robnik-Berry billiard, 
except that the ripples at the bottom of the Mexican 
Hat turn out to be much larger than expected for the 
m = 2 Landau interaction channel. 



II. THE ROBNIK-BERRY BILLIARD 

In this section we will describe how the dot is chosen 
and how the single particle energy levels Eq, and eigen- 
functions (j)a{jc) are determined. We use a trick invented 
by Robnik and BerrvfT^ and elaborated upon by Stone 
and BruusfTsl. Consider a unit circle \z\ = 1 in the com- 
plex plane oi z = x + iy. The analytic function 



w{z) 



z -\- bz"^ + cz^e'-^ 
Vl + 262 3^2 



(7) 



defines a map under which the unit circle in z gets 
mapped into a new shape in w, which will be our dot. 
The shape of the dot can be varied by varying the pa- 
rameters 6, c, and x- The denominator ensures that the 
billiard has the same area (tt) as the unit disc. The wave- 
function (jjaiw, w) = (t>a(u, v) is required to vanish at the 
boundary and obey 
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d d 



4>a{u,v) 



= -4:——4)aiw,w)=ea(t>aiw,w). (8) 

(We have chosen H = 2m = 1). If we now go the 
z plane where the wavefunction is (j)a{w{z),w(z)), the 
Schrodinger equation and boundary condition are 

d d 

-4— — <?!)„(2,2) =£„|w'(z)|20„(z,z) <?!)„(|z = 1|) = 
oz oz 

(9) 

where w' [z) = dw/dz. This differential equation in the 
continuum is next converted to a discrete matrix equation 
by writing 



(t>a{z,z) = 0a(r, ( 



E 



1 



Cf^,(r,0) (10) 



where (r, 9) is the solution to the free Schrodinger 
equation in the unit disk vanishing on the boundary: 



-V2^,(r,0)=72^.,(r, 



(11) 



(That is, these are Bcssel functions in r times angular 
momentum eigenfunctions in 0. ) Feeding this expansion 
into Eqn. Q one obtains the matrix equation 



where 



7» 1] 



(12) 



(13) 



(Without the l/7j in the expansion Eqn. (|10|l . M would 
not have been Hermitian. ) 

In practice one truncates M to a finite size (we used 
585 states) and expects the lower energy levels and wave- 
functions to be unaffected. 
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The parameters b, c, x are chosen to lie in the range 
where classical behavior is chaotic, and where quantum 
chaos as reflected in the eigenvalue distribution has been 
established [Tsl. A value we used repeatedly was b = c = 
.2,x = -85. A nonzero x ensures that the billiard has no 
reflection symmetry. This shape is often called the Africa 
Billiard based on its resemblance to that continent, as 
seen in Fig|21for our chosen values of parameters. 
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FIG. 2: The sfiape of the Robnik-Berry billiard for 6 = c = 
0.2, and x = 0.85. 

We shall refer to these eigenfunctions and eigenvalues 
as exact even though they come from solving a truncated 
problem because we can easily increase the accuracy by 
increasing the size of the truncated Hilbert space. 



g = 2-kKf/{2'k/L) WOF states each of width 2'k/L in 
the tangential direction. Finally Ep = Kp/2m = Kp, 
N ^ kpL^/Ai: and L = VArea = 0r. 

As a test case when we picked the Fermi energy to be 
the 100-th level, we found g = 37. How well is this state 

> at Ep spanned by the g WOF states at the Fermi 
energy? 

First we first take g equally spaced points k„ on the 
Fermi circle and form the WOF states 

i;woF-n{r) = ^e'''--e(dot) (14) 

where 9 (dot) is unity inside the dot and zero outside. 
These states are very close to being orthonormal. For 
example the overlaps of n = 1 state (with k along the y- 
axis) with the others as we go around the circle is shown 
in Fig. ©. 




WOF STATE NUMBER 



A. Testing Assumptions I and II 



Our ability to solve the Schrodinger equation (to high 
accuracy) implies in principle that we can test our two 
assumptions. 

In the next subsection we will test Assumption I, i.e., 
see in detail how well the WOF states serve a basis within 
the Thouless band. 

As for Assumption II, wc and our predecessors [TslIT^ 
have shown that the eigenvalues and single eigenfunc- 
tions obey the distribution expected by RMT for a GOE. 
(The ensemble is generated by varying the parameters in 
w{z)-) ^ 

Similar information about wavefunction correlations is 
not known in the ballistic problem (despite some recent 
progress using supersymmetry methods|l7l|). We did not 
try to do this here since our computing capabilities did 
not allow us to generate an ensemble. 

Instead we computed the fate of the interacting system 
without recourse to Assumptions I and II and compare 
to our predictions based on these assumptions. 

B. Completeness of the WOF basis 



Let Ep be the Fermi energy. Then g ~ y/inN ~ 
\pKEp^ which we arrive at as follows. The Fermi cir- 
cle has a circumference 2'kKp and into this will fit 



FIG. 3: The absolute value of the inner product of 

WOF state number 1 with the otherg — 1 = 36 states. 

Next we ask how much of the state \F > at the 
Fermi energy is contained in the WOF states. We find 
ELi l(" - WOF\F)\^ = .9993. This is a rather re- 
markable result. It says that \F >, which is a vector 
with 580 components (which was the size of our trun- 
cated problem) can be expanded almost completely in 
terms of g = 37 WOF states which are given in advance. 
In other words as one changes the shape of the dot and 
works at fixed Fermi energy, the state \F > changes in a 
random way, but that randomness is only in which par- 
ticular combination of WOF states describes it, not in 
the completeness of the WOF basis. 

While this is very satisfactory we need more to imple- 
ment our scheme: we need to be expand all g states in the 
WOF basis. Here we find that as we move off the center 
of the Thouless band, the fractional norm captured by 
the WOF basis drops. In a typical case, with g = 37, 
there are roughly 12 states (one third of g) where the 
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number lies above 95%. At band edge this drops to 50%, 
as shown in Fig Q). Thus there is inevitably some error 
in transcribing the Landau interaction written in terms 
of the WOF states labelled by k into the basis of g exact 
eigenstates labelled by a. This just means that the loca- 
tion of the critical point will not be correctly predicted 
by our RMT based analysis, as pointed out recently by 
Adam, Brouwer, and SharmafTsj . 



where 




5 10 15 

Exact state index, measured from Ep 



FIG. 4: The norm of the projection of the n exact eigenstate 
from Ef on to the subspace spanned by the WOF states. It 
can be seen that more and more of the exact eigenstate lies 
outside this subspace as one moves away from Ef- 

This concludes our (partial) test of Assumptions I and 
II. We turn to a comparison of our results based on these 
assumptions with a direct solution of the problem with 
no recourse to the assumptions. 



S 



dr 



i/j [idr - Ho 



(18) 



If the fermions are integrated out we will get an effective 
action S^f f (a) . To find the minimum we need just the 
action for static a . In this case it is clear that 



dip dip exp 



drip [idr — Hq + a]ip 



-Eo(a}T 



(19) 

where T — > oo is the length of the imaginary time r- axis 
and Eo{(t) is the ground state energy of ip'^ [Hi:, + cr)ip . To 
find Eq{(j) one simply solves for the single particle levels 
of {Ho+<t) and fills up the ones with negative energy. The 
effective action for static configurations, which is also the 
effective potential, is 



Veff = Eo{a) + 



2u' 



(20) 



At this point we have a mean-field theory. We still 
need to justify its use by showing that fluctuations of 
the collective field cr around its minimum are small. In 
our previous work, based on Assumptions I and II we 
showed that the fiuctuations were indeed small in the 
limit of large since the in front of the actions lim- 
its fluctuations. In the billiard we will justify the mean 
field similarly, based on the depth and curvature of the 
minimum. 

When the Landau interaction is factorized, the hamil- 
tonian whose ground state gives us Eoia) is 



E 

a/3 



(21) 



where, for the case m = 1, for example. 



C. Hartree-Fock solution of interacting problem 

How can the knowledge of the "exact" eigenfunctions 
and eigenvalues in the billiard help in the solution of the 
problem with interactions? The tactic will be illustrated 
in schematic form first. Suppose we have a four- Fermi 
interaction added to a free hamiltonian which in first- 
quantization is given by some differential operator H^. 
Then the path integral becomes 



Z 



dipdipe" 



where 



S= dr 



(15) 



(16) 



Using a Hubbard-Stratanovich transformation we can 
write 



Z = / dipdipdae" 



(17) 



k 



(22) 



and a, /?, k are not restricted to the Thouless band. This 
is because we want to solve the problem without any 
of the assumptions that led to the effective low energy 
theory within the Thouless band. Note that cr has two 
components, because the Landau interaction associated 
with u„i has two parts: 

Vl — (5nk(5nk' (cosTO^k cosm0k'+sinm6'k sinm^k')- 

kk' 

(23) 

Once Seff is known (on a grid of points in the cr plane) 
one can ask if and when the minimum moves off the ori- 
gin. 

So far our considerations have been fairly generic, and 
the Landau interaction has been written in momentum 
space. However, in testing our approach in the billiard, 
we will find it more convenient to represent the Lan- 
dau interaction in real space, since the eigenfunctions are 
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known as linear combinations of Bcssel functions whose 
integrals are best carried out in real space. We have 
carried out calculations for two Laudau parameters cor- 
responding to m = 1 and m = 2. The to = 1 Landau 
interaction is chosen to be (in second-quantized notation) 

x/^^-'*^(-")(5^("*V')(5^*(,^) (24) 

The factors of -j^ — }. on each side of the V have the ef- 
feet of l/|k| in momentum space. Since momentum does 
not commute with the free Hamiltonian iJp, the factors 
have to be placed symmetrically. Note that this corre- 
sponds only to the q ~ Q part of the Landau interaction. 
In reality, all values of q up to the scale El/vf exist in 
the Hamiltonian. Depending on the shape of the dot a 
particular combination of them may break symmetry to 
give the best energy. Still, we expect that since at large g 
we are close to the zero-dimensional limit, the best com- 
bination will consist largely of very small q parts of the 
Landau interaction. In any case, the energy of the true 
symmetry-broken state can only be lower than what we 
calculate, so what we have here is a conservative estimate 
of symmetry-breaking. Similarly the to = 2 interaction 
(also at (f = 0) is 

I Id'r^Hr)j^^;:^{Vl - V^)(5^^*(,-) . 

+-2 I d'r^Hr)j^^^2V,Vy^^^^^{f) 
>< /^^-'*^(-")(2^;ri)T7.2(V0.(V0, 

The integrals are over {w, i«), but can be converted to in- 
tegrals over the disk by using the conformal mapping of 
Eq. lO. Of course, the derivative operators must also be 
transformed in the process. In order to find the matrix el- 
ements of M.af} we had to take the matrix elements of the 
above operators in the basis of exact billiard states. We 
carried out the angular part of the integrals analytically, 
but had to turn to numerical integration to evaluate the 
radial integrals. This is a computationally intensive cal- 
culation, but once the matrix M has been constructed, 
one simply diagonalizes the Hamiltonian of Eq. (|21|l for a 
mesh of cr in the plane, adds up the energies of the lowest 
N particles to obtain the fermionic ground state energy, 
and obtains the effective potential landscape from Eq. 
(|20|l for various values of the coupling strength u. After 
this, it is a simple matter to identify the global minimum, 
which gives us the lowering of ground state energy and 
the value of the order parameter as a function of u. 

Let us proceed to the results, displayed in pictorial 
form. In Fig. we show the absolute value of the 
order parameter, normalized by the nominal value of 
g = \/4nN, for three values of the number of particles 
A'^. The bulk transition happens at u^^n^ = 2. As can be 
seen, there is a nonzero order parameter for any nonzero 



u, and it grows smoothly and continuously as u increases. 
Nothing discontinuous happens at u = 2 or even beyond, 
indicating that the instability does not suddenly become 
first-order at the bulk value of u* . Of course, in these fi- 
nite systems, the Thouless and bulk scales are related by 
a factor g/Air, which is not that large (4.4 for the largest 
system we considered, with A^ = 245). So somewhere be- 
tween u = 2.25 and u = 2.5 the instability seems to reach 
the bulk scale. However, note that the size of systems we 
have considered correspond quite closely to actual ballis- 
tic samples 0, which typically have a few hundred elec- 
trons. Further, the three curves seem to track each other 
fairly closely, indicating that the expectation value of |cr| 
indeed scales with g, as predicted by our earlier work 
based on RMT assumptions. 

6 I , 1 , 1 , 1 , 

N=102l 

5 - - N=245 




u 



FIG. 5: The absolute value of the order parameter normalized 
by (; as a function of coupling strength u for three values of the 
number of particles N. The fact that the curves track each 
other closely indicates that the order parameter does indeed 
scale like g. Furthermore, nothing discontinuous happens at 
the bulk critical coupling strength ul^n. = 2. 

In Fig. we show the corresponding reduction in 
ground state energy normalized by g^. Once again, the 
curves track each other fairly closely, indicating that the 
energy reduction due to interactions is indeed of order 
g^, as predicted by our earlier work. 




FIG. 6: Reduction in ground state energy normalized by 
for three values of the number of particles N. 

In Fig. [7|we show the effective potential landscape for 
m ~ 2, with A^ = 245, at a value of u = 2.15, at which 
the minimum is well-developed, but the order parameter 
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is still within the nominal Thoulcss scale and has not 
reached the bulk scale. The RMT analysis predicted a 
Mexican Hat landscape with "small" ripples (down by 
1/g) in the circle of minima of the Mexican Hat. The 
landscape we see bears no resemblance to this. Instead, 
it appears to be an isolated minimum at a nonzero cr. 
Upon close inspection it can be seen that the minimum 
is shallower in the transverse direction than in the radial 
direction, but this is the only indication we could find of 
a (perhaps) incipient Mexican Hat structure. 
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FIG. 7: Effective potential landscape for symmetry-breaking 
in the m = 2 channel. Instead of a Mexican Hat minimum 
structure with small ripples we see an isolated minimum. The 
minimum does seem shallower in the transverse direction. 

Fig. |S1 shows a similar effective potential landscape for 
symmetry breaking in the m = 1, channel, where the two 
exactly degenerate minima expected from time-reversal 
invariance considerations can be seen. The landscape 
also appears more Mexican-Hat-like than in the m = 2 
case. 



To trace the origin of this difference in behavior, we 
investigated the average absolute value and the 

rms deviation of the matrix element s from the mean 
absolute value, y^(|A/^^p) — (I^'^q/jI)^ for the two cases 

771 = 1,2. The results for the i = 1 (corresponding to 
V^^ for m = 1 and — for m = 2) shown in Fig. |3 
are an energy average for a particular billiard, with the 
parameters b = c = 0.20, J = 0.85. (We have confirmed 
similar behavior of the matrix elements for other param- 
eter values as well.) Fig. shows these quantities as a 
function of the energy difference between the two states 
a and f3. There are two features that are particularly 
noteworthy. 

• There is a "hole" in the m = 1 matrix element near 
zero energy difference. 

• The rms deviation of the m = 2 matrix elements 
from their mean absolute value is huge. As a 
rough estimate, if the matrix elements were Gaus- 
sian distributed complex numbers, the rms devia- 
tion should be roughly half the mean modulus. 
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FIG. 9: A plot of the absolute value and the rms deviation of 
the matrix elements -M^^ from their mean absolute value as 
a function of Sa — ep for the two cases m = 1, 2. 
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FIG. 8: Effective potential landscape for symmetry-breaking 
in the m = 1 channel for N — 200 and u = 2.15. The two ex- 
actly degenerate minima required by time-reversal invariance 
can be seen, as can a Mexican-Hat-like structure. 



However, the i ~ 2 component (corresponding to Vj, 
for TO = 1 and 2Va;Vj, for to = 2) shows very different 
behavior in Fig. ^| While the m = 1 case looks similar 
to the i = 1 component, the fluctuations of the to = 2 
i = 2 component are strongly suppressed by almost an 
order of magnitude below the mean. 

Consider first the "hole" at -Ej? for m = 1. By sym- 
metry considerations alone one can understand that the 
diagonal matrix element Maa for to = 1 has to be zero 
sample by sample in the absence of an external magnetic 
field. Focusing on the x component of the order param- 
eter 

M^a=EcO«(^k)C(k)0„(k) (26) 
k 

By time-reversal invariance ^^(k) — (j)a{—h). Noting 
that 6'_k = fi'k+TT, and that the cos term changes sign, one 
concludes that Maa = —Maa = 0. The reason the "hole" 
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FIG. 10: A plot of the absolute value and the rms deviation 
of the matrix elements from their mean absolute value 

as a function of Sa — £i3 for the two cases m = 1, 2. 

persists for finite energy differences for the operator p = 
— iV can be explained by noting that^^ fo'^ ^ billiard 

p = im[r, H] 

(-iV)a/3 = -im{r)ap{ea - sp) (27) 

which means that the matrix element must vanish at 
least linearly with the energy difference. In fact, such 
"banded" matrix elements have been found for many op- 
erators in ballistic dotspof. 

Consider next the fact that the distribution of the ma- 
trix elements of for the m = 2 case is much broader 
than for the m = 1 case, while the matrix elements 
have a very narrow distribution. The RMT answer would 
have the rms deviation of Map from the mean to be of the 
same order as the mean absolute value. This seems to be 
roughly true for both components of m = 1 but grossly 
untrue for the i = \ component of m = 2. Since it is these 
mesoscopic fluctuations in Maa which determine the size 
of the ripples at the bottom of the Mexican Hat in the 
RMT scenario, this broad distribution of M^p seems to 
be the cause of the failure of the RMT prediction that 
the ripples should be subdominant by 1/g. While it is 
tempting to try to explain this in relation to the shape 
of the billiard (Fig. |2Il, which certainly appears to favor 
an — type of symmetry, a satisfactory explanation 
of the broad distribution of the m = 2, i = 1 matrix 
elements eludes us. 

Our knowledge of the eigenfunctions at the global min- 
imum allow us to compute the effective action for time- 
dependent cr at that minimum. Since the quasiparticles 
couple to this collective field, the interaction induces a 
decay width for the quasiparticles (details can be found 
in ref.Q). In Fig. Illl wc compare the numerically calcu- 
lated values of the width to the parameter-free theoretical 
prediction (solid line) based on RMtQ. On average, the 
RMT based prediction seems consistent with the numer- 
ics, though there is a lot of variation in the widths driven 
by large variations in the matrix elements coupling the 
quasiparticlc levels to the collective mode. 



FIG. 11: A plot of the decay width of quasiparticles induced 
by their coupling to fluctuations of the collective field cr, for 
iV = 245, m = 2, and u = 2.15. The solid line is the theo- 
retical prediction from our previous RMT-based analysis^. 
While the prediction does well on average there are huge vari- 
ations in the widths due to large variations in how strongly 
each level couples to the collective mode. 



III. CONCLUSIONS 



In our earlier work, we used a global RG assumption to 
reduce the problem on the scale of the Thouless energy to 
that of a disordered noninteracting problem with Fermi- 
liquid interactions. This is quite plausible for ballistic 
dots on very general grounds. To proceed further we had 
to make two further assumptions; (i) That the g approx- 
imate momentum states at the Fermi energy were a good 
basis in which to expand the exact disorder eigenstates, 
and (ii) That the wave functions of the exact eigenstates 
in the momentum basis obeyed all the statistical prop- 
erties of RMT. Based on these two assumptions we were 
able to construct a solution to the problem which was 
asymptotically exact in the limit 5 — > oo. This solution 
led to specific predictions for various physical quantities, 
including the size of the order parameter, the reduction 
in energy due to interactions, the shape of the energy 
landscape, and the size of the quasiparticlc decay widths. 

In this paper we have carried out a calculation which 
is still predicated on the validity of the Fcrmi-liquid form 
of the interactions on a scale Ej^ much larger than the 
Thouless energy Et- In retaining this assumption we 
are on firm ground, since after all, the Thouless energy 
can be made as small as one wishes merely by increasing 
the size of the system. We also assumed that the mean- 
field description of the Landau Fermi-liquid interactions 
is valid, which is justified by the fact that the minima 
in the effective potential landscape are indeed of order 
g-5. However, we explicitly eschewed the other two as- 
sumptions that we made in previous work, with a view 
to independently testing their validity. We found that 
our first assumption, that the approximate momentum 
states were a good basis in which to expand the exact 
disorder eigenstates, was extremely good near the Fermi 
energy, but became increasingly inaccurate as one went 
to the edge of the Thouless shell. Wc did not test the sec- 
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ond assumption about wavefunction correlations explic- 
itly, but indirectly through its effects on the predictions 
of our earlier work. We found that most of the predictions 
held up, with the exception of the shape of the effective 
potential landscape in the case of symmetry-breaking in 
the m ~ 2 channel. Even here, the minimum is shaped 
more like a crescent, indicating the possible emergence 
of the Mexican Hat structure at larger values of g (we 
went to the largest value of g that we could given that 
we kept only 585 states and had to keep at least half the 
states empty). We traced the discrepancy back to the 
anomalously broad distribution (compared to estimates 
based on a complex gaussian distribution) of the matrix 
elements M^^. However, we were unable to pin down a 
physical reason for this broad distribution for the m = 2 
case. 

In conclusion, much of the physics we uncovered using 
our RMT assumptions seems to be valid in the Robnik- 
Berry billiard. The second-order transition that we un- 
covered in the f/ — *■ cxd limit seems to indeed be broadened 
into a smooth croossover as expected, belying fears that 
it may be overtaken by a first-order bulk transition. The 
question of how large g has to be before RMT becomes 



fully applicable remains open; another way to phrase the 
question is to ask what the nonuniversal corrections to 
RMT are in ballistic systems. Finally, an important open 
question is whether the broad distributions of matrix el- 
ements of interaction operators is a generic feature of 
ballistic systems, rather than being a special feature of 
the Robnik-Berry billiard, and if so, what physics de- 
termines the width of those distributions. However, our 
results here give us encouragement that the RMT as- 
sumptions can indeed be used with confidence in making 
predictions in ballistic systems, at a qualitative and semi- 
quantitative level. 
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